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Abstract 

Minisuperspace models derived from Kaluza-Klein theories and low energy string 
theory are studied. They are equivalent to one and two minimally coupled scalar 
fields. The general classical and quantum solutions are obtained. Gaussian super- 
position of WKB solutions are constructed. Contrarily to what is usually expected, 
these states are sharply peaked around the classical trajectories only for small val- 
ues of the scale factor. This behaviour is confirmed in the framework of the causal 
interpretation: the Bohmian trajectories of many quantum states are classical for 
small values of the scale factor but present quantum behaviour when the scale fac- 
tor becomes large. A consequence of this fact is that these states present an initial 
singularity. However, there are some particular superpositions of these wave func- 
tions which have Bohmian trajectories without singularities. There are also singular 
Bohmian trajectories with a short period of inflation which grow forever. We could 
not find any non-singular trajectory which grows to the size of our universe. 



PACS number(s): 04.20.Cv., 04.20.Me 



1 Introduction 

One of the main motivations to study quantum cosmology is to investigate if quantum 
gravitational effects can avoid the singularities which are present in classical cosmological 
models If this is indeed the case for the initial singularity, the next step should be 
to find in what conditions the universe recovers its classical behaviour, yielding the large 
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classical expanding universe we live in. In this paper we investigate these problems in 
the framework of minisuperspace models with scalar fields as sources of the gravitational 
field. 

As a first example, we took a non-massive, minimally coupled scalar field, in a 
Friedman-Robertson- Walker universe with spacelike sections with positive constant cur- 
vature. This model can be viewed as an effective multidimensional theory where the scalar 
field is understood as the scale factor of internal dimensions @, or as a Brans-Dicke model 
redefined by a conformal transformation ||. We were able to find the general classical 
solutions. All of them present initial and final singularities. The model is quantized in the 
Dirac way, with arbitrary factor ordering, and the general solution of the corresponding 
Wheeler-DeWitt equation is found. To interpret the solutions, we first adopted the 'peak 
interpretation', where a prediction is made when the wave function is sharply peaked in a 
region and almost zero outside this region || . A gaussian superposition of WKB solutions 
was constructed. By employing the stationary phase condition, we were able to show that 
this superposition is sharply peaked around the classical trajectory only for small values 
of the scale factor. Hence, contrarily to what is usually expected, the classical limit is 
recovered for small values of the scale factor, not for large ones. A consequence of this fact 
is that the initial classical singularities continue to be present at the quantum level. In 
order to confirm this strange behaviour, we also adopted an alternative interpretation of 
quantum mechanics which was not constructed for cosmology but which can be easily ap- 
plied to a single system: it is the causal or the Bohm-de Broglie interpretation of quantum 
mechanics [[|. It is completely different from the others because it is an ontological in- 
terpretation of quantum mechanics. In the case of non-relativistic particles, the quantum 
particles follows a real trajectory, independently of any observations, and it is accompa- 
nied by a wave function. The quantum effects are brought about by a quantum potential, 
which can be derived from the Schrodinger equation. It is a rather simple interpretation 
which can be easily applied to minisuperspace models ||. In this case, the Schrodinger 
equation is replaced by the Wheeler-DeWitt equation, and the quantum trajectories are 
the time evolutions of the metric and field variables, which obey a Hamilton- Jacobi equa- 
tion with an extra quantum potential term. The application of this interpretation to 
some of the quantum solutions of our problem shows exactly the same behaviour as found 
previously: the Bohmian trajectories behave classically for small values of the scale factor 
while the quantum behaviour appears when the scale factor becomes large. Singularities 
are still present. However, when we make superpositions of these wave functions, the 
initial singularity disappears for some special cases, but none of these special trajectories 
grows to the size of our universe. 

The other case studied involves two minimally coupled scalar fields. They can be 
viewed as a tree level effective action of string theory where the second scalar field comes 
from the Kalb-Rammond field [0j. They can also be understood as generalized Brans- 
Dicke type models, which can be derived from compactification of multidimensional theo- 
ries with external gauge fields || . The results obtained in this case were analogous to the 
preceding one. Along the lines of the peak interpretation, gaussian WKB superpositions 
predicts a classical universe for small values of the scale factor because they are peaked 
around the classical trajectories in this region. Adopting the causal interpretation to in- 
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vestigate the singularity problem, we found, as before, that many of the solutions present 
classical behaviour when the scale factor is small (and hence singularities) but behaves 
quantum mechanically when the scale factor becomes large. 

This paper is organized as follows: in the next section we describe the classical min- 
isuperspace models of both one and two scalar fields models, presenting their general 
classical solutions. In section 3 we quantize these models obtaining their corresponding 
Wheeler-DeWitt equations and their respective general solutions. In section 4, the gaus- 
sian superpositions of WKB solutions are constructed and their peak along the classical 
trajectories are exhibited. In section 5, the causal interpretation of quantum cosmology 
is shortly reviewed and applied to the quantum solutions. We end with comments and 
conclusions. 



2 The Classical Models 

Models with two scalar fields that interact non-trivially between themselves can be ob- 
tained from different theoretical contexts. Considering Kaluza-Klein supergravity theo- 
ries, keeping just the bosonic sector, and reducing to four dimensions, leads to effective 
actions with gravity plus two scalar fields, one of them coupled non-minimally to the 
Einstein-Hilbert lagrangian; the two scalar fields have an interaction between them. More 
generally, every time we consider multidimensional models with gauge fields, and reduce 
them to four dimensions, we find such structure. String theories, in particular, have an 
effective action in four dimensions given by the expression, 

L = ^e-^R + <P w r J -^H, uX H^ , (1) 

where is a dilaton field and H^ u \ is a Kalb-Ramond field which in four dimensions is 
equivalent to a scalar field £ 

In order to keep contact with this variety of models, all of them having great impor- 
tance in high energy conditions, we will consider the general lagrangian 



L = y/=g\^>R-u^--^-) , (2) 

where u is a coupling constant. We remark the non-trivial interaction between and 
£. For the string effective action, uj = — 1 and for Kaluza-Klein theories u = where 
d is the dimension of internal compact spacelike dimensions. If we perform a conformal 
transformation such that g^ = 0~ 1 <? AU ,, we obtain the lagrangian 
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(3) 



where the bars have been suppressed. From Eq. (0) we deduce the field equations, 
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+ = o , (5) 

(p K(p 

n ^_ 2 iipp = , (6) 



3 
2' 

We consider now the Robertson- Walker metric 



where k — u> + ~. 



r/.s' J - -y'dt 2 + -^ 1 [dr 2 + r{ t W 2 - ,nr((-))d./)}. (7) 
1 + 4 r 



where the spatial curvature e takes the values 0, 1,-1. The equations of motion are, for 
N = l, 



,d 2 3e ^/^\2 _|_ jjvj^ 
V a 2 2> j 2>' 



+ 3^0-^ + ^- = , (9) 
a (p K(p 

a ■ d) ■ , 

£ + 3 _£_ 2 r£ = o . (10) 

a <p 

We will be interested in the case e = 1. In what follows we will consider separately the 
cases where £ = const, and £ ^ const.. 



2.1 One scalar field minimaly coupled to gravity 

Henceforth, we consider in Eqs. (|8|j9| JT0|) £ = constant. The solutions of the resulting 
equations can be easily found if we reparametrize the time coordinate as dt = a 3 d6. The 
integration procedure is standard, and we just give the final results: 




In these expressions, A, B and C are integration constants. The universe expands from 
an initial singularity untill a maximum size and then contract to a final singularity. Note 
that a oc ts for small a. For A = 1 and B = C, we obtain the implicit relation: 



2.2 Two scalar fields coupled to gravity 



Considering the fields (j> an d £ in equations (|| and using again the same parameter 

9 as defined previously, we find the following solutions: 

£ = ^ + £ Ktanh (M) ; (14) 
* ~ ^^#) 1 <15) 

5 1 (16) 



cosh(y§|C|(0 + E)) 
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In these expressions A, B, C, D and E are constants. The qualitative behaviour of the 
scale factor is the same as in the preceding case (compare fll2|) with (|16D). Again we have 
a oc ts when a is small. For A = 0, B = C,D = E and k = |, we can find a simple 
implicit relation between a, <fi and £: 

0(a) = p^6a 2 ; (17) 

«(0 = V^l^-f)^ ; ( 18 ) 

m = ^-\e ■ (19) 



These implicit classical relations, together with Eq. fll3|) , will be compared with the 
trajectory on which the semi-classical wave function of the corresponding quantum model 
is peaked. 

3 Quantum Solutions in Minisuperspace 

We return to the lagrangian (g) and we insert on it the metric (^). The action takes the 
form 

S = [ Ldt (20) 



where 

12ad 2 .„ „ .a?(b 2 „a 3 £ 2 „„„ r 

L =— - {3+2 ^w- 2 w-~ UNa ■ (21) 

From ( PH) we obtain the conjugate momenta, 

ad 

*a = 24— , (22) 



^ = -2(3 + 2^)0 , (23) 



« = ~ A w ■ (24) 



We can now construct the hamiltonian H, which takes the form 

r 2 



H = N 



7T„ 



J.2„2 



48a 4(3 + 2uj)a 3 8a 3 



V 2 

f + 12a 



NH . 



(25) 



Variation of N yields the first class constraint 7i ~ 0. The Dirac quantization procedure 
yields the Wheeler-DeWitt equation by imposing the condition: 



and performing the substitutions 



7T„ 



7T t 



ws> = 



<9 2 p 9 



<9a 2 a <9a 
(9 2 g 9 



d 2 



(26) 

(27) 
(28) 
(29) 



where p and g are ordering factors. We have set h — 1. The Wheeler-DeWitt equation in 
the minisuperspace reads 



a 
12 



2 r 



^ aa ~ * a 

a 



(3 + 2a>) 



2 



(30) 



where V$(a) = 48a 4 . 

We will solve this equation for the cases £ = (one scalar field) and £ 7^ (two scalar 
fields). 



3.1 Solutions with one scalar field 

Discarding the field £, we have to solve the equation, 



a 2 




2 




12 


(3 + 2a;) 


a 





(a) * . 



(31) 



Supposing the separability of this equation, we can write \l/(a, 0) = a(a)/3((f>) leading to 
two ordinary differential equations for a and (5: 



p 

«aa + - = V a (a)a , 
a 



where 



y„ = 12 ( 48a 2 - 



W) = -(3 + 2a;)- 



'02 



(32) 
(33) 

(34) 



a 




Figure 1: Behaviour of a(a) for n G I for the one scalar field case, with p — 1, k — 1, 
A a = 1 and B a = 0. The dashed and continous lines represent the imaginary and real 
parts of a respectively. 



k being an integration constant. The solutions for a and (5 are, 



a k (a) = a^' 2 



AJ n {\2a 2 ) + B a K n (12a 2 ) 



(35) 



i (l+m- ? )/2 

^(0) = A^ 1 -^)/ 2 + 5/ , (36) 



m 



with n = V(p |) 4 ^ an( j m _ ^(g — 1)2 _ 4(3 _|_ 2cj)A;. The function a does not exhibit 
an oscillatory behaviour unless n e I. For this case, a oscillates for small values of a, 
increasing or decreasing for large values of a, suggesting that a classical phase may occur 
for small values of a only. The function (3 has an oscillatory behaviour, for all 0, if m e I. 
In figure 1 we show the behaviour of the real and imaginary parts of a for n G I. The 
complete solution of the Wheeler-DeWitt equation is 

V(a,<l>)=jA(k)a k (a)l3 k (<l>)dk . (37) 
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3.2 Solutions with two scalar fields 



For the case where both scalar fields are non null, the Wheeler-DeWitt equation in the 
minisuperspace reads, 



a 

12 



* aa ~ * a 

a 



(3 + 2cj) 

We use again the separation of variables method writing \l/(a, 0, £) 
tion separates in two: 



Kp(a)* . (38) 
E(a,0)A(£). Equa- 



A 



-8fciA 



a 
12 



1 — 'ao 1 — a 

a 



(3 + 2u) 



V 3 (a,<f>)E 



(39) 
(40) 



with Vs(a,4>) = 48a 4 — 4/c!0 2 , fc x being an integration constant. Writting S(a, 
a(a)/3(0), we obtain two ordinary equations: 



£ 

a 



with 



V^(a)a , 

14(a) = 12(W-^ 

(3 + 2u) Uk x - % 



The solutions for a, (3 and A are 

Q-M 2 \AJ n {l2a 2 ) + B a K n (12a 2 



(41) 
(42) 

(43) 
(44) 



a\a) 



n 



a 



(p - l) 2 - A8k 2 



a-ff)/a 



V ro ^(3 + 2^)^ ^3 + 2^x0^ 



ni 
A(0 



(g- l) 2 -4(3 + 2w)fc 2 



(45) 



(46) 
(47) 



The coefficients A's and S's are constants. The general solution of the Wheeler-DeWitt 
equation is 

V(a,<f),£) = J A{k 1 ,k 2 )a k2 (a)(] kuk2 (<f))*k 1 (0dhdk 2 . (48) 

In general, a is an exponentially growing or decreasing function of a. If the order of the 
modified Bessel functions is imaginary, a may exhibit an oscillatory behaviour. However, 
for these cases, a oscillates for small values of a, increasing or decreasing for large values 
of a, suggesting again that a classical phase may occur only for small a. This behaviour 
is displayed in figure 2. 
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a 




a 



Figure 2: Behaviour of a(a) for the two scalar fields case for p — 1, &2 = 1, A a = 1 and 
-B a = 0. The real part is represented by the continous line while the imaginary part is 
represented by the dashed line. 
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4 The WKB Approximation 

One way to try to obtain the transition to the classical regime from the quantum solutions 
is to employ the WKB approximation, like in usual quantum mechanics. This is achieved 
by rewriting the wave function as, 

tt = exp(^) , (49) 

substituting it into the Wheeler-DeWitt equation, and performing an expansion in orders 
of h in S, 

s = s + hs 1 + h 2 s 2 + ... (50) 

The classical solution must be recovered by constructing a wave packet from So: 

* = J A(k )exp(^-S )dk , (51) 

where ko is an integration constant. As in the preceding sections, we will analyse the 
WKB approximation separately for the cases with one and two scalar fields, respectively. 



4.1 WKB approximation with one scalar field 

In this case, we have S = S(a, 0), and the WKB expansion in the minisuperspace Wheeler- 
DeWitt equation, leads to the following equations connecting S and Si. 



<r_ //AV 2 
12 \~da~ 



2 fdS G 



2 r 



a 
12 



d 2 S \ fdS ^ fdS 
da 2 ) 
d 2 S 



3 + 2u\d(p 

ip ( dSo 



+V*(a) = 



V da J V da 
dS \/dSi 



-2 



a V da 
^iq_fdS 



3 + 2u 

First we get a solution for S . It can be obtained by taking, 

S (a,<f>) = S (a) + S (<P) , 
leading to two differential equations: 

'dS (a^ 2 



da 



12['^-48o 2 

^ a 



(52) 

(53) 

(54) 

(55) 
(56) 



where ko is a separation constant. These equations admit the following analytic solutions: 

^3(k — 48a 4 ) — ^/3/c arctanh( 



S (a) 



lk -48a 4 \ 



k 



-) 



+A 



±^/(3 + 2oj)koln(f) + Bo , 
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(57) 
(58) 



where An and Bn are integration constants. We follow the same procedure in order to 
obtain a solution for Sx(a,<p), considering first Si(a,0) = Si (a) + Si(0). We get the 
solutions, 



Si (a) 
Si(0) 




arctanh 



'kn -48a 4 



+i 



P 



— 1 i 

— In a + - ln(48a 4 - jfe ) + A 1 



2u 



kn 



In. 



+B 1 



,(59) 
(60) 



where Ai and B\ are integration constants. From the solution for So (a), we can easily 
see that only for kn > we can obtain an oscillatory behaviour of the wavefunction for 
small values of a, while for k < the wavefunction has an exponential behaviour for any 
value of a. Similarly, if (3 + 2u)k n > 0, then exp [^So(0)] is oscillatory for any value of 0, 
otherwise it has an exponential behaviour. Hence, for kn > and u > — |, exp ^S o (a,0) 
oscillates for small values of a and any value of 0. 

We can construct a wavepacket from the above solutions through the expression, 



A{k Q ) exp [-S (k , a, <f))]dk 



(61) 



where the function A(kn) is a sharply peaked gaussian centered in fc > with width a. 
Examining Eq. (|57D, we can see that So (a) becomes very large when a becomes very 
small. Hence, in the integral (|6lD , constructive interference happens only if 

dSn(a, <f> 



dkn 







(62) 



which implies a relation between kn, a and 0, kn = kn(a, <fi). The wave function turns out 
to be: 

*(a, 0) = A[k (a, 0)] exp{^S [A;o(a, 0), a, 0]} . (63) 

As the gaussian is sharply peaked at ko(a, 0) = ko, then we obtain that the wave function 
(HQ is sharply peaked at k (a, 0) = k . It can be verified that this relation is exactly the 
classical relation (0) with k playing the role of the integration constant A. 



4.2 WKB approximation with two coupled scalar fields 

We follow the same procedure as before, writing the wave function \1/ in terms of S(a, 
and performing an expansion in orders of fi. The final equations for So and Si are: 

a 2 fdSnV 2 /<9S \ 2 2 /<9Sn 



12 V da 
12 

3 + 2cj 



3 + 2uj V <90 
<9 2 So\ /dS ( 



da 2 



d 2 Sc 



da 
dS 



2 

dSi 
da 

V d6 



d£ 



•Vnf(a) 







-0, 



(64) 



| ip fdS 
a V da 

( 



^iqfdSn 



+ 



^2 r 



■9 2 Sp 

de 



dS \/dS 1 



(65) 
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Imposing again the ansatz So (a, <fi, £) 
equations: 

<95 n (a)\ 2 



5*0(0) + 5*0(0) + 5o(£); we obtain the following 



da 

fdS (<f>)Y 



dS (Q 



(3 + 2u) 
2k , 



48a' 
K 



- h 



(66) 
(67) 
(68) 



where K Q and k$ are separation constants. The solutions are: 
5 (a) 



So(<f>) 
So(0 



± 



± 



3(K - 48a 4 ) - V3fsT arctanh 



iKn - 48a 4 



K 



+A 



(3 + 2lu)(K - k (f) 2 ) - V(3 + 2cu) J fs:oarctanh 



-Bn 



±j2k £ + C , 



where Aq, B$ and Co are integration constants. As in the one scalar field case, we can 
find solutions for Si but they are not important for the construction of the wave packet 
in our approximation. The solutions 5o will be enough to recover the classical trajectory. 
First we note that K > leads to a oscillatory behaviour for exp [^5o(a)]. On the other 
hand, if (3 + 2u)K > 0, keeping K > 0, then exp [|5 (</>)] is oscillatory for any value 
of <p when k < 0, or only for small values of <fi when k Q > 0. If (3 + 2uj)K < 0, then 
exp [^5o(0)] has an exponential behaviour for any value of <fi when k$ < or for small 
values of <fi when ko > 0. 

We consider now the superposition given by 



*(a, 0, f ) = J J A(k , K ) exp ^5 (a, 0, f , k , K )dk dK 



(69) 



where A(ko, Kq) is a bidimensional gaussian function, centered on fco > and Kq > with 
width o\ and cr 2 , respectively. As before, 5 (a) becomes very large for small a. Hence, 
we have to guarantee constructive interference by the condition, 



/ <95o(a,0 ,Q 
V dk 



\k =k 



+ 



<95o(a,0,Q 

dK lKo=Ko 







(70) 



The implicit relations coming from ([T(]) are the same as the classical relations (p~T| , p~8| , [T9| ) . 
The classical limit is again recovered only for small a. 



5 The Perspective of the Causal Interpretation 

In this section, we will apply the rules of the causal interpretation to the wave functions 
we have obtained in section 3. We first summarize these rules for the case of homogeneous 
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minisuperspace models. In the case of homogeneous models, the supermomentum con- 
straint 7i l is identically zero, and the shift function iVj can be set to zero without loosing 
any of the Einstein's equations. The hamiltonian is reduced to general minisuperspace 
form: 

H GR = N(t)H(p a (t),q a (t)), (71) 

where p a (t) and q a (t) represent the homogeneous degrees of freedom coming from IP- 7 (a;, t) 
and hij(x,t). The minisuperspace Wheeler-De Witt equation is: 

H(p a (t),q a (t)Wq) = 0. (72) 

Writing \1/ = Rexp(iS/h), and substituting it into (f72|), we obtain the following equation: 



where 



1 r)1 81 

2 Mq » ) WaW? +U{q " ) + Q{q " ) = 0, (73) 



Q(q») = —fapj^-, (74) 



R dq a dqf3 



and fapiq^i) and U(q fl ) are the minisuperspace particularizations of the DeWitt metric 
Gijki and of the scalar curvature density — h 1 ' 2 R( 3 > (hy) of the spacelike hyper surf aces, 
respectively. The causal interpretation, applied to quantum cosmology, states that the 
trajectories q a {t) are real, independently of any observations. Eq. (|73|) is the Hamilton- 
Jacobi equation for them, which is the classical one ammended with a quantum potential 
term ([F4]), responsible for the quantum effects. This suggests to define: 

P a = (75) 
dq a 

where the momenta are related to the velocities in the usual way: 

f = s+™SL . (76) 

To obtain the quantum trajectories we have to solve the following system of first order 
differential equations: 

dS(g a ) = fa(3 ]_dqp (77] 
dq a N dt 1 ' 

Eqs. (f77|) are invariant under time reparametrization. Hence, even at the quantum level, 
different choices of N(t) yield the same spacetime geometry for a given non-classical 
solution q a (t). There is no problem of time in the causal interpretation of minisuperspace 
quantum cosmology. Let us then apply this interpretation to our minisuperspace models 
and choose the gauge N — 1. 



13 



5.1 One scalar field 

The general solution of the Wheeler-DeWitt equation is given by 

tf(a,0) = j A{k)a k {a)f3 k {4>)dk . 

where 



i-p 

a fc = a 2 



A Q / n (12a 2 ) + 5 a K n (12a 2 



^ < , 1 - m -1 Br l + m-q 

p k = Arf—s- + -?-4>-^- 

m 



with 
and 

The momenta are 



n = 



(p - l) 2 - 48k 



m = ^J(q - l) 2 - 4(3 + 2u)k 



ir a = 24aa , 

tt^ = -2(3 + 2co)a 3 -^ . 

The causal interpretation states that the momenta are also given by 

dS{a,<f)) 



7I> 



da 
dS(a, 0) 



(78) 

(79) 
(80) 

(81) 
(82) 

(83) 
(84) 

(85) 
(86) 



where S(a, <fi) is the total phase of the wave function Hence, the Bohmian trajectories 
will be solutions of the following system of equations: 



24ad = 



dS(a, 



-2(3 + 2u)aT 



da 
dS(a, 0) 



2 



d<t> 



(87) 
(88) 



The quantum potential for this problem can be calculated in the usual way. We substitute 
^ = Re lS into the Wheeler-DeWitt equation, obtaining the Hamilton- Jacobi like equation 
with the extra quantum potential term Q: 



a^(dS_\ 



4 



12\daJ Z + 2u\d<t>J 



f^) 2 -48a 4 + Q = , 



where 



Q = 



R 



a 2 fd 2 R pdR\ 



fd 2 R q dR\ 



12\da 2+ ada) 3 + 2cA<90 2 + d<j> ) 



(89) 



(90) 
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Let us apply this interpretation to the simplest case $ = a>k{a>)/3k{4>)- Then the wave- 
function has the form, 

* = J R 1 (a) J R 2 (0)e i[5l(a)+52W] , (91) 

since S(a,<fi) = Si (a) + S 2 (0). This implies that flSTD becomes independent of 0. From 
Eq. (|9~0|), we see that Q(a,0) = Qi(a) + Q2(4>)- To simplify the calculations we set 
Ap = Bp = 0, and p = q = 1. We will first calculate the dynamics of the scale factor 
when a is small, in order to see if there are singularities. In this approximation we can 
take just the first term of the series represantion of I n (x) |l( 



00 1 / T \n+2l 



f^l\T(n + l + l) V2 



For n real, the modified Bessel function I n (x) is real and the phase of is zero. Hence, the 
Bohmian equation flS^) yields that a is a constant. It is a nonsingular quantum solution 
but with little physical interest. Hence, in a first moment, we will restrict ourselves to 
the case where n is a pure imaginary number. Combinations of these two situations will 
be analyzed afterwards. 

In the case where n is pure imaginary, can be written as 



a k = Co x iu = c e iulnx , n = iu = ±iV3k . (93) 
The phase and the norm are, 

Si(a) = v\nx , ( . 



Defining x = 12a , Eq. fl87|) becomes 

dS v 



whose solution is 



dx x ^ ^ 



ut.i , 
—\* . (96) 



v 4 ' 

We can see that, in accordance to what was suggested in previous sections, the behaviour 
of the quantum trajectory is like the classical one for small a, and the singularity will still 
be present. Note that if v is positive we have expansion, while if v is negative we have 
contraction. For the scalar field, we have: 



k = A p <tr = Bj^ , (97) 



where u = — y (3 + 2uj)k. The phase and the norm are 

S 2 (<j>) = uln0 , 



# 2 (0) = B a . ' 
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From Eq. fl8~8]) we have, 



Using Eq. (RBT) we get 



3 + 2w 3 - OS u 

2 — a *=oj = -^ ' (99) 



(/) = tV3m^) 5 (ioo) 

which is also the classical behaviour for A = 1. 

It is not surprising that we have obtained the classical behaviour. Since R\(a) and 
i?2(o) are constants, the quantum potential is zero, and there is no quantum effect. Note 
also that solutions (pBj) and ( |100| ) satisfy the hamiltonian constraint with V(a) = 48a 4 
neglected because a is very small. 



For very large a, the Bessel function I n (z) can be approximated to [10 

e 12a 2 

4(x)~ Cl , (101) 

a 

where C\ is a complex constant. In this case, we have: 

5i(a) = const. , (102) 

12a 2 

RJa) = led . (103) 

a 

The quantum trajectory is evidently 

a = const. , (104) 
which is not the classical one. For the scalar field, Eqs. ( p8| , |98| , p!04D now yields, 



2 V 3+2ui 



(105) 



which is the classical behaviour for <p in this regime (note that as a = const., 6 oc t in 
equation (|Tl~D). The different behaviours of the scale factor and the scalar field can be 
explained with the quantum potential. For the scale factor, Qi will be given (see Eq. (|90|) 
with p = 1): 

Qx(a)=48a 4 + ^ (106) 

which is of the same size of the classical potential V(a) = —48a 4 , and hence responsible 
for this quantum behaviour. For the scalar field, as R2{<P) is constant, Q2{4 > ) = and 
the scalar field continues to follow its classical trajectory. Note that the hamiltonian 
constraint is also satisfied in this limit. Hence we have again obtained the strange result 
where the classical limit happens only for small values of a. 

Let us now take some superpositions of the afc(a) and Pk{4>) given in Eqs. (|79| , |80D . 
For definitiness, we will choose p = q = 1, and u = 0. We will continue to take only pure 
imaginary n's. Combinations of real and pure imaginary n's do not change qualitatively 
the results. The wave function will be given by, 

* = 4, (12a 2 ) U<T U > + Bi<p\ (107) 

i=l L 
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Figure 3: Field plot of a versus using the causal interpretation for the superpositions 
of the wave functions in the one scalar field case, in the region of small a. 
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Figure 4: Plot of a particular singular solution for a(t), coming from figure 3, which begins 
with an inflationary phase. 
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where 



A, 
Bi 



n ! 



1 

3 

b - e 



t.3 

A 2 
B 2 



1 




, Mi = 

« 2 = n 3 
A 3 = 
B 3 = 1 



Using Eq. (|92"D we obtain, for small a: 



12a' 



'12a 2 



u 2 = u 3 



(108) 



12a 2 



(109) 



Figure 3 shows a field plot of a versus 0, for small a, for the Bohmian equations (|S 
with S being the phase of the wave function (|109| ). We can see that there are periodic 
solutions with very small oscillations around a. They are eternal quantum universes 



which never grows. The solutions which grow beyond the validity of ( |109|) are singular 
and inflationary. This can be seen in figure 4. 

This result suggests that the initial singularity can be avoided only if we superpose 
eigenfunctions of opposite frequencies. However, it seems to be be difficult to obtain, in 
scalar field models, non-singular universes with long expansion period [O]. 



5.2 Two Scalar Fields 



In this case, we have studied the quantum trajectories driven by wave functions obtained 
from Eq. ( [48"D for some particular A{k\, k 2 ). 



i) We have fixed p = q = 1, uj = 0, k\ = —jn, k = k 2 , A, 



A, 



0, B n 



B n = A> 



B\ = 1 and A(ki, k 2 ) = §<5(^i + To) tanh(7r^/3^2). Using a result of Ref. HlOl , we obtain 



= cosh(^y|e)^^tanh(7ry3l)^(x)^ 



V3k' 



icf))dk 



arctan(— ) 

x 



(110) 



The quantum trajectories can be calculated from the following equations (in the gauge 
N=l): 



The solutions are: 



tt„ — 24ad — 9 d s a — x 2 + ^2 , 

Qa 3 <p _ dS_ _ x 2 +<j> 2 +x 



1 a 
7T</> 



ACK _ dS _ Q 
^ 4> 2 ~ d£ ~ ' 



x z +4> 2 



[III) 



In 



In 



„, , -62- 

£ = const. , 



(112) 
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where 77 = / — is the conformal time and C is an integration constant. For small a, when 
r] approaches ± 2 ~ , these functions tends to: 



a(t) oc ts ; 
4>(t) oc oc a 2 
£ = const. 



(113) 



which is exactly the classical behaviour for uj = 0. When a is not small, the trajectories 
are not classical (compare (|112| ) with Eqs. (|TJ,|I^,|n|)). This can be seen by inspecting 
the quantum potential. For two scalar fields it is given by 



^ ~ RlV2 Ufa? + a~da)~3 + 2uj 



d 2 R q dR 
^7 + -^ 



d 2 R 



2 oe 



For our particular problem, we obtain 

1 (x 4 - 20 2 x + 4 ) 



Q 



x 2 + 



(114) 



(115) 



For small a we have Q oc a 2 (remember that <fi oc a 2 in this limit). In this domain, the 

2 '<9S\ 2 1 x 2 (j) 2 



kinetic terms dominate: 



a 2 n 2 



12 



a 

12 V da J 

k2 



3(x 2 + 



m2 



x 



oc const. 



<P 2 



X 



X 2 + 



oc const. 



(116) 
(117) 



Hence the quantum potential becomes negligible when compared with the classical 
kinetic terms. For a not small, for instance, when n — > yielding a — > a max and — > 00, 
the quantum potential diverges, 

Qoc0 2 , (118) 



while the classical potential and kinetic terms behave like 



oc 



oc 



i> 2 

l2 



Vd oc a 4 



(119) 

(120) 
(121) 



Hence, together with K^, the quantum potential becomes the more important term. This 
behaviour of the quantum potential explains why the trajectories are classical for small a 
and quantum otherwise. 



ii) Let us now take p = q = 1, u 
B a = Bp = A x = 1 and A(k t , k 2 ) 
the following wave function (see Ref. [|Kj): 

7r 2 f x , , 1 



0, &i = i k = k 2 , A a = At 



0. 



5{k\ — I) sinh(7r 1/3^2)^ /3fc-(-y/2e J 4 ). Then we obtain 



exp 



x 

5' 



+ 



- & 



cxp 



3 x 



:i22) 
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The equations of motion are: 



VT a = 


24ad = 


- dS - 
da 


= 24^# 


7T0 = 


° p 


as 

90 


_ _j_ 

X 


7T ? = 


* 02 


as 

<9€ 


[2 

V 3 



(123) 



The solutions are, 



1 

a = £3 



= Coa 2 = c t3 , (124) 
C = -|ci|ts + C 2 • 



Note again that these solutons approach the classical one when a is small but are com- 
pletely different when a is large. Once again, this can be explained by the behaviour of 
the quantum potential when compared with the kinetic and classical potential terms. The 
quantum potential is given by 

6 2 

Q = 4&a i -\ . (125) 
The kinetic and classical potential terms are given by, 

lv a 3x 2 ' 
K, - -J£- 

7_0^ 2 ' (126) 
V a = -48a 4 . 

For small a, Q, V a and goes to zero while K a and are constant. For large a, 
Q, and K% become comparable and large, while K a and continue to be constant. 
Hence, in this situation, the quantum potential becomes important, driving the quantum 
behaviour of the Bohmian trajectories. Note that the sum of ( |125| ) with (|126|) gives zero 



because the Bohmian trajectories must satisfy the hamiltonian constraint ammended with 
the quantum potential term. We have also calculated the Bohmian trajectories for other 
exact wave solutions of the Wheeler-DeWitt equation. All of them present the same 
behaviour. 

We conclude this section by stating that in quantum cosmology it is not necessary 
that the classical behaviour appears when a is large, while quantum behaviour is present 



when a is small. It can indeed be the reverse. This was already pointed out in [12] and 
we presented specific examples illustrating this fact. It should also be commented that 
the result of this seciton using the causal interpretation are in qualitative agreement with 
the results of the previous section. 



6 Conclusion 

In this paper, we have studied classical and quantum minisuperspace models containing 
one and two scalar fields. We have shown that all classical solutions are singular. After 
quantizing the models, we have obtained the general solution of the Wheeler-De Witt 
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equation. Usually, the solutions are oscillatory when the scale factor is small and not os- 
cillatory when the scale factor becomes larger. This suggests that non-classical behaviour 
may occur when the sacale factor is large. We studied gaussian superpositions of WKB 
wave functions to investigate if they correspond to quasi-classical states, as suggested 
in Ref. [|13j. We have shown that indeed these wave functions are peaked around the 
classical trajectories in configuration space, but only for small a. 

After, we applied the causal interpretation of quantum mechanics to these models. In 
this interpretation, it is possible to calculate quantum trajectories, independently of any 
observations. We have shown that the trajectories calculated following this interpretation 
usually present the classical behaviour when the scale factor is small and non-classical 
behaviour when a is large, as suspected. This means that these quantum trajectories 
still presents an initial singularity. We have also seen that if we superpose negative with 
positive frequency solutions, then we can find trajectories which are no more classical 
for small a. We can have eternal periodic quantum universes with very small oscilations. 
These universes, however, never scape the Planck length. There are also singular solutions 
with a short period of inflation which grow forever. We could not find any non-singular 
solution which grows to the size of our universe, with a classical limit for large a. 

The fact that quantum behaviour happens when a is large is not surprising. It was 
already obtained in Ref. [II] and suggested to exist when the scale factor grows like t 1 ' 3 
in Ref. |12|, which is our case. Hence, in quantum cosmology it is not necessarily true 
that large scale factors implies classical behaviour. For the scalar field models we have 
analyzed in this paper, the reverse seems to be more usual. It means that it is possible to 
have in our universe some degrees of freedom which still behave quantum mechanically in 
spite of it being very big. This gives us some hope of being possible to detect or experience 
quantum cosmological effects in the real universe we live in, bringing quantum cosmology 
to the realm of testable physical theories. The problem should be to find which degrees 
of freedom can possess this property. To do this, we should improve this minisuperspace 
model with the accretion of small inhomogeneous perturbations, which contain an infinity 
number of degrees of freedom, and see what happens with the new inhomogeneous degrees 
of freedom. We will get closer to the real universe but we will have to face new technical 
and interpretational problems. This will be the subject of our future investigations. 
Acknowledgements We thank CNPq (Brazil) and CAPES (Brazil) for financial support. 
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